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We show that a recent non Born-Oppenheimer calculation of dipole moments exhibits obscure 
points and is not consistent with the well known Hellmann-Feynman theorem. 
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CN . I. INTRODUCTION 

X> 

|"T i . Some time ago Cafiero and Adamowicz[l] calculated the dipole moments of LiH and LiD without resorting to 

\& • the Born-Oppenheimer (BO) approximation. Their approach consists in the expansion of the eigenfunction of the 

Coulomb Hamiltonian in a basis set of floating s-type explicitly correlated Gaussian functions 1]. The dipole moments 

' t— | '' obtained in this way are essentially identical to the experimental values [2] and the authors claim that their calculations 

simulate experiment more closely than any previous onesQ. They also applied basically the same approach to the 

^3 | calculation of molecular polarizabilities^, 

The purpose of this paper is to analyze that non Born-Oppenheimer calculation of dipole moments of diatomic 

'— 1 . molecules in the light of the Hellmann-Feynman theorem[6j . In Section|TT]we outline some well known properties about 

the molecular Hamiltonian that are necessary for the remaining sections, in Section IIIII we develop the Hellmann- 

' Feynman theorem for optimized variational wave functions, in Section IIVI we apply it to the calculation of dipole 

^ | 1 moments described by Cafiero and Adamowiczfl] and in Section [V] we draw conclusions and summarize the main 

CO | results of this paper. 
00 ' 

o 

oo 
O 



II. PRELIMINARIES 



In this section we outline some well known properties of the molecular Hamiltonian that are necessary for present 
discussion of the nonadiabatic calculation of dipole moments. A more detailed analysis was provided earlier by other 
authors [?1 Q]. Following Cafiero and Adamowicz [lj we consider a nonrelativistic Coulomb Hamiltonian of the form^O] 



^ 2m, ^ ^ 47reor i7 - V ' 

i i j>i J 

where rrii, qi andpi are the mass, charge, and momentum, respectively, of particle i located at position r», = |r, — r 3 | 
is the distance between particles i and j, and eo is the vacuum permittivity. The first step in the treatment of such 
system of particles is the separation of the motion of the center of mass which we carry out by means of the linear 
transformation 
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so that the Hamiltonian operator JTJ) becomes 

— 1 +EE 7^ = ^ + Hcm (3) 

where p[ is the momentum conjugate to rj. In this way we remove the motion of the center of mass given by Hqm and 
are left with the translation-free molecular Hamiltonian Hm provided we choose the transformation matrix elements 
tij conveniently 

Following Cafiero and Adamowicz[l| we can choose = ycm and r'j = rj — r±,j > 1, where ycm and ri denote the 
positions of the center of mass and the heaviest nucleus in the molecule, respectively. Other authors have proposed 
a more symmetric transformation: = Tcm and r' = rj — r/vcM> where y^cm is the position of the center of mass 



of the nuclei [9|. In this case the coordinates of the nuclei are not independent and one of them can be expressed in 

n n 

terms of the others [S|. One advantage of the transformation chosen by Cafiero and Adamowicz[l] is that the form of 
the potential-energy function remains simple because Vji = |r^| and ry = |r^ — r£|, i,j > 1. 

The inversion operator i changes the sign of the coordinates and momenta of all the particles: iTji = —Vj, 
ipjl" 1 = —f)j, where i _1 — I. It is clear that the translation-free Hamiltonian operator is invariant under inversion 
iHm^ 1 = Hm- Therefore, if ip and E are an eigenfunction and its eigenvalue, respectively, of Hm then iHm^ — 
Huiip = Eiip. If the state is nondegenerate we conclude that iip — ±ip because i 2 — 1 (the identity operator). In 
particular, the ground state is nondegenerate and spherically symmetric 0, [kJ. 

The classical dipole moment for a neutral distribution of charges is given by the well-known expression 



j 

It follows from i^mT 1 = —fx and what was said above that (ip\fj,\ip) = for any nondegenerate molecular state. In other 
words, the outcome of the quantum-mechanical nonadiabatic calculation of the dipole moment as the expectation 
value of the corresponding operator is zero for any molecule in its ground state. In order to circumvent this difficulty, 
Cafiero and Adamowicz considered the classical interaction between the dipole moment and the field e thus obtaining 
the semiclassical Hamiltonian 

H e =H M -e-Li (5) 

that behaves as iff e (e)£ -1 = H e (~e). If we choose e along the z axis, e = (0, 0, e), then e • /i —t[i z . If ip e and E t are 
an eigenfunction and its eigenvalue, respectively, of H e we have iH e (e)ij; e = H e (—e)iip e = E e iip e . Consequently, if the 
state is nondegenerate then the energy is an even function of the field: E e (—e) = E e (e). 

Throughout this paper we omit the fact that the Hamiltonian operator ((5|) does not support bound states because 
it was not an issue in the nonadiabatic calculation of molecular dipole moments and polarizabilities[ll 3, 



III. THE HELLMANN-FEYNMAN THEOREM 

In what follows we briefly discuss the Hellmann-Feynman theorem for variational wave functions [6| because the 
nonadiabatic calculation of dipole moments has necessarily been based on them. Consider a trial function $ that 
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depends on a set of adjustable parameters a = {ai, a^, ■ ■ ■ , a„}. If we differentiate the variational energy E(a) given 
by 



E($\ $) = <$|£/"|$) 



with respect to the variational parameters and require that 



7T> = »■'='■-' 

oat 



then the optimal trial function satisfies 

|i H-E\^)^ opt + (^\H-E 



d<f> 
da; 



= 0,i = l,2,. 



(6) 
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Suppose that the Hamiltonian operator H depends on a parameter A (which may be a particle charge or mass, the 
intensity of an applied field, etc). The optimized variational function will depend on A through the optimal parameters 
a opt . Therefore, if we differentiate equation ([6]) with respect to A 



dE 
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H - E \<f>) + (<f>\ H - E 



(9) 



and take into account the variational expression then we prove that the optimized variational function satisfies 
the well known Hellmann— Feynman theorem [f]]: 



dE _ dH 
d\ ~ \ ~d\ 



(10) 



where (^Aj — ($|^4|<f>) /($|<f>). Of course, any eigenfunction of H and its corresponding eigenvalue satisfy this 
relationship. 



IV. NONADIABATIC CALCULATION OF DIPOLE MOMENTS 

In order to calculate the dipolc moment of the diatomic molecule Cafiero and Adamowicz 1] fitted three points of 
the energy curve E e (e) with a second-degree polynomial and obtained \i z from the coefficient of the linear term. 

E e (e) «e + eie + e 2 e 2 + ... (11) 



They resorted to polynomials of higher degree in other calculations 
According to the Hellmann-Feynman theorem (| 10(1 we expect that 

9 -§ = - (12) 

for the kind of variational function that Cafiero and Adamowicz [l] chose to solve the Schrodinger equation approxi- 
mately. Consequently, 

dE 

^7 = ei = - <A.)„ = (13) 

unless the variational function is not fully optimized or it does not recover the correct symmetry as e — > 0. This result 
is consistent with the argument of Section |H] that E e (e) is an even function of e. 



4 



The simple arguments outlined above clearly show that a nonadiabatic calculation, either as an expectation value 
or by fitting the energy with a polynomial function of the field intensity should produce the only result of zero dipole 
moment. However, Cafiero and Adamowicz 1] managed to obtain the dipole moments of the diatomic molecules 
by means of the following "trick" "The energy was calculated for each basis with three electric field strengths, 
e z = 0, —0.0016, and —0.0032 a.u., and the energy curve was fitted with a second order polynomial in ej. /i is then 
the first order coefficient of this fit." In fact, the theoretical dipole moments estimated in this curious way appear to 
converge smoothly towards the experimental ones as the number of basis functions increases [l]]. Later they explained 
this procedure in a more detailed way[4|: "We use the finite field approach in the present work, i.e., we calculate 
the energies of the system for several field strengths, we fit the energy as a function of the field strength with a 
polynomial, and last we use the polynomial to determine the energy derivatives with respect to the field at the zero 
field. As is clear from the discussion above, the non-BO energy of a molecule at the field strength / is identical to 
the energy at the field — /, because when the direction of the field changes the orientation of the molecule follows the 
field direction. Thus, for any system in the H2 isotopomer series the energy is an even function of the field, and if it 
is approximated by a polynomial, only even powers need to be used. This obviously results in a zero dipole for any 
system if the dipole moment is determined as described above. An alternative approach is to apply stronger fields and 
only use energies calculated for positive field strengths in generating the polynomial fit. In this case the use of both 
even and odd powers is appropriate. As we have shown in our previous work on LiH[l|, the dipole moment derived 
from our non-BO calculations with the procedure that uses only positive fields and polynomial fits with both even 
and odd powers match very well the experimental results. Thus in the present work we will show results obtained 



using interpolations with even and odd-power polynomials." Summarizing: Cafiero and Adamowicz [lj, |4| tells us that 
it is possible to fit an even function by means of a polynomial one with even and odd powers of the variable provided 
one chooses sufficiently great values on the positive side of the variable axis. It is clear that the nonzero coefficients 
of the odd powers should be the product of numerical errors. Otherwise, if we apply this argument to, for example, 
cos(/), then we will prove that sin(0) 7^ 0. A question arises: if the outcome of the linear term in the expansion 
of the energy as a function of the electric-field intensity was due to numerical errors, why did the dipole moments 
estimated by Cafiero and Adamowicz [l| agreed so accurately with the experimental ones?. Notice that this same 
"trick" led to negligible dipole moments for the A2 diatomsjj] which is consistent with the classical view that they 
should have zero dipole moments: "Applying the same approach to homonuclear species (H2, D2, and T2) should give 
the dipoles identically equal to zero. In our calculations, these actually come out to 10 -8 . This small noise which 
entered our calculations was due in part to the previously mentioned fact that the zero-field wave function we use is 
not an eigenfunction of J 2 as it should be. The level of noise introduced is negligible, as 10 -8 is 4 orders of magnitude 
smaller than the size of the dipole moments for the heteronuclear species." 

It seems that Cafiero and Adamowicz 1] chose suitable trial wave functions for the calculation of the dipole moments 
that did not reflect the expected symmetry of the exact eigenfunctions at zero field strength (otherwise the linear term 
of the polynomial (fTT|) would have been zero) . The success of their calculation was probably due to the particular 
placement of the floating Gaussian functions that reflects the ionic character of the chemist's classical picture of the 
molecules: [l| "Thus the centers corresponding to the hydrogen nucleus were scattered from about 2.9 to about 3.1 
bohrs. The lithium nucleus was, of course, placed at the origin of the internal coordinate system. The functional 
centers corresponding to the electrons were located primarily on the two nuclei, with two electrons at the origin (about 



5 



0.0 ± 0.001 bohrs in all three directions) and two electrons near the H nucleus (about 3.05 ± 0.06 bohrs) per basis set. 
This reflects the strong ionic character in the lithium/hydrogen bond. The LiD non-BO wave function was optimized 
starting from the converged LiH wave function." Therefore, it is not clear that the variational wavefunction reflects 
the correct symmetry of the system for the values of the electric-field intensity chosen for the fit. If the variational 
functions were allowed to approach the symmetry of the exact eigenfunctions of H e ^o = Hm as the field vanishes, the 
linear terms of the fitting polynomials would be zero as well as the predicted dipole moments 10( ■ This may probably 
be the case of the homonuclear diatoms Q|. 

V. CONCLUSIONS 

When the field is on, the wavefunction no longer have spherical symmetry, the expectation value of /t is nonzero 
and depends on the field strength e. However, it is well-known that E e (—e) — E e (e) and the expansion of the energy 
about e = should exhibit only even powers of the field strength (as we have briefly shown in Section [TTJ) . In other 
words, the linear term should be zero in agreement with the Hellmann-Feynman theorem (| 13[) . 

We have shown that if one allows the optimized variational wave function to recover the correct spherical symmetry 
when the field intensity vanishes then the linear term of the polynomial (fTTj) should be zero as well as thepredicted 
dipole moment. Apparently, the optimized variational wave function used by Cafiero and Adamowicz [l[ 0, Q 
retains, when e — > 0, the cylindrical symmetry it has for e > 0. This may be due to the particular placement of the 
floating gaussians in space. Another hint was given by Cafiero et al[ll|: "Since the basis functions used here are 
eigenfunctions of J z but not J 2 , they are of correct symmetry only for the non-zero-field calculations. In order to 
maintain equal representation of all the points in the fit, though, they are used for the zero field point as well. This 
should introduce only mild contamination by excited rotational states and has not been a practical problem in the 
past." However, that contamination may be the cause of the occurrence of the linear term in the fit of EL (e). For that 
reason, we believe that the apparent success of the approach proposed by Cafiero and Adamowicz [l[ [jj [j, Q should 
be further investigated. 

It is clear from the following referee's comment that this problem is not well understood: "There is no spherical 
symmetry when the electric field is included, thus there is no inconsistency in the calculations of Cafiero and Adamow- 
icz in this respect." However, from the arguments already given above it is clear that the coefficient of the linear 
term should be the expectation value of jl z at zero field strength that obviously corresponds to a spherical-symmetric 
ground-state. The referee is one of the supporters that sin(0) ^ 0. 
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